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Multiwavelength studies of radio relics at merger shocks set powerful constraints on the relics origin and formation mech- 
anism. However, for X-ray observations, a main difficulty is represented by the low X-ray surface brightness far out in 
the cluster outskirts, where relics are typically found. Here, we present XMM-Newton results from a 130-ks observation 
of CIZA J2242. 8+5301, a cluster at z = 0.19 that hosts a double radio relic. We focus on the well-defined northern relic. 
There is a difference of ~ 55% between the temperature we measure behind the relic, and the temperature measured with 
Suzaku. We analyse the reasons for this large discrepancy, and discuss the possibility of reliably measuring the temperature 
beyond the northern relic. 
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1 Introduction 

Radio relics are diffuse, Mpc-scaled, arc-shaped sources, 
typically found at large cluster-centric radii. They are 
strongly polarized, with magnetic fields roughly aligned 
with the length of the relic (for reviews, see e.g., Briiggen 
et al. 2012, Feretti et al. 2012). The commonly-accepted 
theory explaining the origin of relics is that they trace 
cluster- merger shocks. The shocks will accelerate electrons 
and ions to relativistic energies, producing a population of 
synchrotron-emitting particles observable at radio frequen- 
cies. The shocks will also heat and compress the ICM ther- 
mal plasma, and these effects should be observable at X-ray 
wavelengths in the form of temperature and density jumps. 
The Mach number of a shock can be calculated from the 
spectral index at the front of the relic, assuming a certain ac- 
celeration mechanism and pre-existing particle population. 
At the same time, it can be calculated from the X-ray prop- 
erties of the gas on both sides of a relic, using the Rankine- 
Hugoniot conditions. Therefore, combined radio and X-ray 
observations of relic-hosting galaxy cluster mergers allow 
us to shed light on the origin of radio relics by testing corre- 
lations between the thermal and non-thermal effects of clus- 
ter merger shocks. 

So far, only two merger shocks have been X-ray- 
confirmed at radio relics (Abell 754, Macario et al. 2011; 
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Abell 3667, Finoguenov et al. 2010), even though the num- 
ber of known radio relics is much larger (~ 50, Feretti 
2012). The major difficulty of X-ray observations is that 
relics are often detected at large distances from the cluster 
centre (Vazza et al. 2012), where the ICM surface brightness 
is very low. Even if cluster emission is present at these dis- 
tances, modelling the spectra is complicated by confusion 
with the background. Therefore, for correctly measuring the 
thermal properties of the gas and the shock Mach number, it 
is essential to accurately model the sky background, and to 
characterize the systematic errors introduced by uncertain- 
ties in the instrumental and sky backgrounds. 

The sky background is composed of thermal foreground 
Galactic emission, which dominates at soft energies, and 
of the cosmic X-ray background (CXB; emission from un- 
resolved sources), which dominates at high energies. The 
CXB is well-described by a power-law with a photon index 
of 1.41 (De Luca & Molendi 2004). At high Galactic lati- 
tudes, |6| > 20°, the foreground is typically modelled as the 
sum of Local Hot Bubble emission (unabsorbed, Tx ~ 0.1 
keV) and Galactic Halo emission (absorbed, Tx ~ 0.2 
keV). Close to the Galactic plane, the hydrogen column 
density increases, and the foreground emission becomes 
more complex. Indeed, additional thermal or line compo- 
nents often need to be included in the foreground model. 
Masui et al. (2009) detected additional foreground emission 
from a Suzaku observation in the Galactic direction (255°, 
0°), which is best-described by a thermal component with 
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T ~ 0.75 keV. Simionescu et al. (2011), in their analysis 
of Suzaku observations of the Perseus cluster (b w 15°), 
found a thermal foreground component with T pa 0.6 keV. 
George et al. (2009) measured densities and temperatures in 
PKS 0745-191 beyond the virial radius (r2oo; the radius at 
which the mean enclosed density is 200 times the critical 
density of the Universe at the redshift of the cluster), but 
their results were challenged by Eckert et al. (2011), who 
discovered a discrepancy between the Suzaku and ROSAT 
surface brightness profiles, most likely attributed to an ad- 
ditional foreground component that was not taken into ac- 
count by George et al. (2009). Later Suzaku observations 
near PKS 0745-191 allowed for an improved background 
analysis (Walker et al. 2012), showing the need for an ad- 
ditional thermal foreground component with T k, 0.6 keV; 
interestingly, the background fit was not improved by the 
inclusion of the Galactic Halo component. These exam- 
ples serve as cautionary tales about the importance of ac- 
curate background modelling, in particular when analysing 
the faint outskirts of galaxy clusters. 

At the redshift of CIZA J2242.8+5301, 1 arcmin corre- 
sponds to 192 kpc. All the errors are given at the la level. 

2 CIZA J2242.8+5301 

CIZA J2242. 8+5301, a merging cluster at z = 0.19, hosts 
a spectacular northern radio relic: 2-Mpc long, ~ 75-kpc 
wide, strongly polarized (> 50%), and with magnetic fields 
aligned with the relic (van Weeren et al. 2010). The merger 
takes place in the plane of the sky (van Weeren et al. 201 1), 
which reduces the number of unknowns that can be intro- 
duced by projection effects. Therefore, CIZA J2242. 8+5301 
is as a textbook example for studying the connection be- 
tween radio relics and merger shocks. 

We have successfully proposed for XMM-Newton and 
Chandra observations of CIZA J2242. 8+5301 . The main re- 
sults of these observations (130 and 200 ks, respectively) 
will be discussed in Ogrean et al. (2012) and Ogrean 
et al. (in preparation). CIZA J2242. 8+5301 was also ob- 
served with Suzaku (Akamatsu & Kawahara 2011; here- 
after, AK1 1). Worryingly, there appear to be significant dif- 
ferences between the XMM-Newton- and Suzaku-measmed 
temperatures, which further affect the X-ray-derived shock 
Mach number at the northern relic. Here, we present the 
XMM-Newton spectral analysis of the sky background in 
the direction of the cluster and the ICM properties across 
the northern relic, and discuss possible reasons for the dis- 
crepancy between the XMM-Newton and Suzaku results. 

3 Spectral analysis 

We reduced the data using the XMM-Newton Extended 
Source Analysis Software (ESAS), integrated in the Scien- 
tific Analysis System (SAS) vl2.0.1. The details of the data 
reduction can be found in Ogrean et al. (2012). 




Fig. 1 XMM-Newton MOS+pn 0.5 - 4 keV exposure- 
corrected image, with the instrumental background sub- 
tracted. Overlaid in green are the regions used for extract- 
ing spectra. The dashed moon-shaped region shows the de- 
fined background. In cyan are WSRT 1.4 GHz radio con- 
tours, drawn at [1, 5, 10, ...] x 100 ^Jy/beam. 

Spectra were extracted from the regions shown in Fig- 
ure[H which are all located within 12 arcmin of the detector 
centre. The background region is beyond r2oo ~ 2 Mpc 
(Ogrean et al. 2012) from the merger centre, and should 
therefore contain only minimal cluster emission. For each 
region, we also created instrumental background spectra. 
Source and sky background spectra were grouped to a min- 
imum of 30 counts per bin. 

Modelling the background in the direction of the cluster 
is complicated by the cluster's low Galactic latitude (b ~ 
—6°). The best fit to the background is given by a model 
consisting of CXB, LHB, GH, and hot foreground (HF; Ma- 
sui et al. 2009, Simionescu et al. 2012) emission (see Sec- 
tion 0). We used APEC models for the thermal components 
(LHB, GH, HF), and a POW model for the CXB compo- 
nent. Absorbed components (GH, HF) were described by 
WABS*APEC models, with coupled, but free, hydrogen col- 
umn densities. The GH temperature was fixed to 0.08 keV, 
the average of the values found by Sidher et al. (1996) and 
Kuntz & Snowden (2000). The hydrogen column density 
was frozen to 3.34 x 10 21 cm -2 , corresponding to the 21-cm 
line measurement of the Galactic hydrogen column density 
(Dickey & Lockman 1990) at the position of the cluster. 

Cluster emission in all regions was modelled as a single- 
temperature, absorbed thermal component (WABS*APEC). 
The temperature and normalizations were free in the fit, but 
coupled for the MOS and pn spectra extracted from identical 
sky regions. The metallicity was fixed to 0.2 solar, and the 
cluster redshift to 0.1921 (Kocevski et al. 2007). 

To maximize the number of constraints, the sky back- 
ground and target spectra were fitted simultaneously. The 
instrumental background was substracted before fitting. The 
fits were performed in Xspec vl2.7.1, using Anders & 
Grevesse (1989) abundances, and Balucinska-Church & 
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Table 1 Best-fit background parameters. Unstarred values were obtained using a background model that includes HF 
emission. Starred values were obtained with a background model without HF emission, as used by AK1 1 (see Section|4]l. 
Temperatures are in units of keV. Power-law normalizations are in units of photons keV -1 cm~ 2 s _1 arcmin -2 at 1 keV. 
The normalizations of the APEC components are in units of cm~ 5 arcmin -2 . Values marked with a dagger (f) were fixed 
in the fit. 
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Table 2 Best-fit temperatures for the partial annuli in the 
northern sector, with (Tx) and without (T^) a HF compo- 
nent included in the background model (see Section|Ui. 
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McCammon (1992) photoelectric absorption cross-sections. 
The fit had a X 2 of 1985.38, with 1833 degrees of freedom. 
The best-fit parameters are summarized in Tables Q] and [2] 
The temperature beyond the relic is not constrained, hence 
we do not include it in Table [2] 



4 Discussion and conclusions 

Unfortunately, the XMM-Newton data does not allow us to 
constrain the temperature beyond the northern relic. The ra- 
dio spectral index is —0.6 ±0.05, which in the linear regime 
corresponds to a (putative) shock Mach number of 4.6+q'q. 
In the X-ray, we do not see a temperature or density jump 
across the relic. However, assuming the relic traces a shock 
of Mach number 4.6, a post-shock temperature of ~ 15 keV 
implies a pre-shock temperature of ~ 2 keV. AK11, using 
Suzaku observations, measured a pre-shock temperature of 
1.7 ± 0.20 keV, although they do not detect a surface bright- 
ness jump across the relic; additionally, they measured a 
post-shock temperature of 6.7 ± 0.45 keV, more than two 
times lower than the XMM-Newton temperature measured in 
approximately the same region. In fact, the highest temper- 
ature in the Suzaku temperature profile across the northern 
relic is ~ 8 keV. One major difference between our analy- 
sis and that of AK1 1 is the background modelling. Unlike 
us, they did not include a HF component. We modelled the 
XMM-Newton spectra using the same model as AK11, to 
see if the different background model can account for the 
temperature discrepancies. For this, we simply needed to re- 
move the HF component from the model described in Sec- 
tion [3] The results of the sky background and ICM spectral 
fits are summarized in Tables Q] and [2] 



An F-test was run to check whether it is truly reason- 
able to add a HF component to the background model. This 
yields a probability of ~ 10~ n . The low value of the prob- 
ability indicates that, based on the XMM-Newton data, the 
background is more accurately described by a model that 
includes a HF component. 

When using the same background components as AK1 1, 
the XMM-Newton temperatures are significantly lowered, 
especially at large distances from the merger's centre, and 
we also derive a pre-shock temperature that is consistent 
with that of AK1 1. However, the XMM-Newton and Suzaku 
results are still not consistent in the post-shock region. 
Another candidate for the temperature discrepancy could 
be residual soft proton (SP) contamination of the XMM- 
Newton event files. We checked the SP-filtered event files 
for residual SP contamination using the FlN_OVER_FOUT 
routine (De Luca & Molendi 2004), in the energy bands 
8-12 and 10 - 12 keV for MOS, and 12 - 14 keV for 
pn event files. The event files were diagnosed as not con- 
taminated by SP. 

The difference between the XMM-Newton and Suzaku 
temperatures could be blamed on calibration uncertainties. 
We note that at the time of the Suzaku observation and when 
the paper of AK1 1 was published, there were some uncer- 
tainties in the calibration of XIS0 and XIS 1 detectors, intro- 
duced by the rising contamination layer thickness of XIS0, 
and the change in the charge injection level of XIS1. Un- 
fortunately, we are not able to tell if these uncertainties are 
enough to explain the temperature discrepancy. 

Another problem in measuring temperatures far out in 
the cluster outskirts is the systematic errors introduced by 
the background model. This could constitute an issue for 
measuring the pre-shock temperature at the northern relic. 
Leccardi & Molendi (2008) analysed a sample of ~ 50 
hot (T > 3.5 keV) galaxy clusters at intermediate red- 
shifts (0.1 < z < 0.3) from the XMM-Newton archive, and 
showed that for source-to-background count rate ratios be- 
low 0.6 in the energy band 0.7 — 10.0 keV, the temperature 
can be biased downwards by a factor of a few. In Figured 
we show the temperature profile across the northern sector 
as a function of source-to-background count rate ratio in the 
energy band 0.5 — 7.0 keV. For the background, we assumed 
the background region shown in Figure [T] The source-to- 
background count rate ratio is less than 0.5 in the pre-shock 
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Fig. 2 Top: Temperature and normalization profiles 
across the northern sector in Figure Q] using a background 
model with and without HF emission (squares and cir- 
cles, respectively). Bottom: Temperature profiles plotted 
against the source-to-sky background (left) and source-to- 
total background (right) count rate ratio. The profiles for 
the background model including HF emission are slightly 
shifted, for clarity. 

region. Compared to the threshold ratio set by Leccardi & 
Molendi (2008), our ratio is probably even less encouraging, 
given that the threshold in the 0.5 — 7.0 keV band should be 
higher than in 0.7 — 10.0 keV (due to lower background). 
Therefore, it is not surprising that we are not able to prop- 
erly constrain the pre-shock temperature from our XMM- 
Newton observation. At the moment, no study similar to that 
of Leccardi & Molendi (2008) was done for Suzaku, so it is 
difficult to say how much the temperature measurement in 
the pre-shock region is affected by background systematics. 
In conclusion: 

1 . When modelling the sky background at low Galactic 
latitudes, it is essential to consider all possible foreground 
components. In the case of CIZA J2242.8+5301, if the HF 
component is not included in the spectral model, the post- 
shock temperature is biased downwards by ~ 30%. If the 
northern relic traces a merger shock, then the measurement 
of the colder pre-shock temperature will be less affected by 
the incorrect background model, as the HF and ICM temper- 
atures are more similar in the pre-shock region; hence, this 
will also bias downwards the X-ray-derived Mach number. 

2. Using the presented XMM-Newton observation, it is 
impossible to measure the temperature beyond the north- 
ern relic in CIZA J2242. 8+5301 while keeping systematic 
uncertainties under control. This might be possible with 
Suzaku, due to its lower instrumental background. How- 
ever, lacking statistical studies about the necessary source- 
to-background count rate threshold for Suzaku spectral anal- 
yses, it is currently hard to tell if the Suzaku pre-shock tem- 
perature would be constrained or significantly biased down- 
wards even when using a correct background model. 

3. Neither XMM-Newton, nor Suzaku detect a surface 
brightness discontinuity across the northern relic. In the 
case of XMM-Newton, this is due to the high background 
level, with the ICM emission beyond the relic being consis- 
tent with zero. However, this does not appear to be the prob- 
lem in the Suzaku brightness profile. So far, all the large- 



scale relics with X-ray-confirmed shocks show a surface 
brightness jump across the relic. While not necessarily im- 
possible (substructure along the line of sight could, in prin- 
ciple, complicate the brightness profile and smooth density 
jumps), it would be puzzling if the northern relic in CIZA 
J2242. 8+5301 was an exception. From the paper of AK11, 
it is not clear whether the Suzaku surface brightness profile 
was background-subtracted. If it was, the lack of a bright- 
ness discontinuity could indicate, for example, small-scale 
CXB fluctuations; if, let's say, the CXB is lower in the post- 
shock region than the assumed average, this would cause 
the post-shock surface brightness to be lower than that de- 
rived, and it would also bias the ICM temperatures towards 
smaller values. 
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